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Abstract. Radio timing observations of miUisecond pulsars (MSPs) in support of Fermi LAT 
observations of the gamma-ray sky enhance the sensitivity of high-energy pulsation searches. With 
contemporaneous ephemerides we have detected gamma-ray pulsations from PSR B1937H-21, the 
first MSP ever discovered, and B1957H-20, the first known black-widow system. The two MSPs 
share a number of properties: they are energetic and distant compared to other gamma-ray MSPs, 
and both of them exhibit aligned radio and gamma-ray emission peaks, indicating co-located 
emission regions in the outer magnetosphere of the pulsars. However, radio observations are also 
crucial for revealing MSPs in Fermi unassociated sources. In a search for radio pulsations at the 
position of such unassociated sources, the Nanfay Radio Telescope discovered two MSPs, PSRs 
J2017H-0603 and J2302H-4442, increasing the sample of known Galactic disk MSPs. Subsequent 
radio timing observations led to the detection of gamma-ray pulsations from these two MSPs as well. 
We describe multiwavelength timing and spectral analysis of these four pulsars, and the modeling 
of their gamma-ray light curves in the context of theoretical models. 

Keywords: gamma rays: observations, pulsars: millisecond, pulsars: individual (B1937H-21, 
B1957-I-20, J2017-H0603, J2302-H4442) 
PACS: 95,97 

INTRODUCTION 

In contrast with normal pulsars, most millisecond pulsars (MSPs) have binary orbits, 
making them extremely difficult to find in blind searches of the gamma-ray data alone, 
because of the vast parameter space to be searched and the sparsity of the gamma-ray 
photons. Gamma-ray pulsation searches from MSPs therefore rely heavily on observa- 
tions at radio frequencies, with observations lasting only a few minutes to several hours. 
For known MSPs, searches for pulsations are done by folding the high-energy data using 
astrometric, rotational and orbital parameters measured through radio timing observa- 
tions. Fourteen MSPs have been observed to emit pulsed gamma rays with this strategy 
in the first two years of data recorded by the Fermi Large Area Telescope (LAT) [1, 2, 3]. 
Another approach consists in searching Fermi unassociated sources for radio pulsations 



from unknown pulsars, a strategy which has so far led to more than 20 discoveries of 
MSPs at the Green Bank, Parkes, Nangay and Effelsberg telescopes [3]. 

The famous MSPs J1939+2134 (a.k.a. B1937+21) and J1959+2048 (B1957+20), re- 
spectively the first millisecond pulsar ever discovered [4] and the first known black- 
widow system [5], are examples of MSPs observed to emit pulsed gamma rays by fold- 
ing the data using timing parameters measured at the Nangay and Westerbork telescopes 
[6]. The Nan9ay Radio Telescope has also contributed to the wealth of discoveries of 
new MSPs in Fermi unassociated sources by finding the two previously unknown MSPs 
J2017+0603 and J2302+4442 [7]. 

TWO MSPS WITH ALIGNED RADIO AND GAMMA-RAY 

EMISSION 

All MSPs so far detected with the Fermi LAT in gamma rays have spin-down energy 
loss rates E above 10^^ erg s ^ making J1939+2134 (E = 1.1 x 10^^ erg s"^) and 
J1959-I-2048 (E = 7.5 x 10-^^ erg s^^) likely gamma-ray emitters. Nevertheless, the two 
MSPs are located at low Galactic latitude, and hence suffer from intense contamination 
by the diffuse gamma-ray emission in the Galaxy. In addition, both have relatively 
large distances of d = (7.7 ± 3.8) kpc for J1939-I-2134 [8] and d = (2.5 ± 1.0) kpc 
for J 1959-1-2048 [9], so that their "spin-down luminosities" E/d^ are on the low end 
of gamma-ray MSPs detected thus far with Fermi. 

Both MSPs have been monitored at radio wavelengths in the context of the timing 
campaign of best gamma-ray candidates for detection by Fermi [10]. We phase-folded 
the first 18 months of gamma-ray data recorded by the LAT using pulsar ephemerides 
based on radio timing data taken at the Nangay and Westerbork radio telescopes, and de- 
tected pulsed gamma-ray emission from J 1939-1-2 134 and J 1959-1-2048. Figure 1 shows 
radio and gamma-ray light curves of J1939-I-2134. The close alignment of radio and 
gamma-ray emission peaks, also observed for PSR J 1959-1-2048 and J0034— 0534 [2], 
suggests that for this subset of gamma-ray MSPs low and high-energy radiation is 
produced in the same regions of the magnetosphere. As illustrated in Figure 2, radio 
and gamma-ray profiles of PSRs J1939-I-2134 and J1959-I-2048 are well reproduced by 
"altitude-limited" Two-Pole Caustic (TPC) and Outer Gap (OG) models, in which both 
radio and gamma-ray emission peaks result from caustics generated at high altitude with 
respect to the light cylinder [11], confirming that radio and gamma-ray emissions seem 
co-located for these MSPs. 

NEW MSPS DISCOVERED AT NAN^AY 

Of the 1451 gamma-ray sources listed in the Fermi Large Area Telescope First Cata- 
log (IFGL) [12], 630 are not associated with counterparts known at other wavelengths. 
Sources of gamma-ray emission with no known counterparts and with emission proper- 
ties that resemble those of gamma-ray pulsars (lack of flux variation over time, spectral 
cutoff at 1 - 10 GeV) can be searched for pulsations from radio pulsars. Observations of 
Fermi unassociated sources likely to be pulsars conducted at Nangay yielded the discov- 




FIGURE 1. Radio and gamma-ray light curves of PSRs J 1939+2 134 and J2302+4442. Two pulsar 
rotations are shown for clarity. See [6, 7] for details on the analysis of the radio and gamma-ray data. 
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FIGURE 2. Model fits for the radio and gamma-ray light curves of PSR J 1939+2 134. 



ery of two MSPs in binary orbits, PSRs J2017-I-0603 (P = 2.896 ms) and J2302-I-4442 
{P = 5.192 ms). Because the two MSPs were found in unassociated gamma-ray sources, 
they were likely to emit gamma rays in a pulsed way. We verified this assumption by 
folding the gamma-ray data using ephemerides based on timing data collected at the 
Nan9ay, Jodrell Bank, and GBT telescopes, and detected pulsations in the gamma-ray 
data recorded after the discoveries. We then used the entire Fermi dataset to extend the 
timing solutions to the past, thereby improving the measurement of astrometric and ro- 
tational parameters. Detailed timing parameters for the two new MSPs are given in [7]. 

Figure 1 shows an integrated radio profile at 1.4 GHz and gamma-ray light curves 
of PSR J2302-I-4442. These light curves are reminiscent of those observed for other 
gamma-ray MSPs, with relatively complex radio profiles and two gamma-ray peaks 



offset from the radio emission. Detailed modeling of the observed radio and gamma-ray 
light curves using standard outer magnetospheric models shows that, similarly to other 
gamma-ray MSPs, the gamma-ray emission seems to originate at high altitudes in their 
magnetospheres (see [13] for details on the light curve modeling of PSR J2017+0603). 

CONCLUSION 

Most of the >20 MSPs newly discovered in unassociated gamma-ray sources are ex- 
pected to be pulsed gamma-ray emitters, so that more than a third of currently known 
gamma-ray pulsars are MSPs, making them a prominent class of Galactic sources of 
gamma rays. The four MSPs discussed here illustrate the importance of radio observa- 
tions in support of the Fermi mission, such as pulsar timing measurements and pulsar 
searches conducted at the Nangay radio telescope. The Nanfay observatory will con- 
tinue to perform timing observations of MSPs for Fermi and search for new pulsars in 
unassociated sources. A third MSP, PSR J2043-I-1711 (P = 2.380 ms) has been discov- 
ered recently in a Fermi source with no counterpart, and will hopefully be followed by 
others. 
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